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QCD rescattering mechanism for diffractive deep inelastic scattering 

Roman Pasechnikp Rikard EnbergjH and Gunnar Ingelmaro 
Department of Physics and Astronomy, Uppsala University, Box 516, SE-751 20 Uppsala, Sweden 

We present a QCD-based model where rescattering between final state partons in deep inelastic 
scattering leads to events with large rapidity gaps and a leading proton. In the framework of this 
model the amplitude for multiple gluon exchanges is calculated in the eikonal approximation to all 
orders in perturbation theory. Both large and small invariant mass Mx limits are considered. The 
model successfully describes the precise HERA data on the diffractive deep inelastic cross section in 
the whole available kinematical range and gives new insight into the density of gluons at very small 
momentum fractions in the proton. 



I. INTRODUCTION 

Hadronic processes with a hard scale involved consti- 
tute an indispensable tool for probing the QCD dynamics 
of quarks and gluons, and through the QCD factorization 
theorems [lj that separate physics at small and large dis- 
tances, one may also study the dynamics of soft processes 
with small momentum transfers. Hard quark and gluon 
interactions at small distances are thus not affected by 
soft interactions and are described in perturbative QCD. 
The most problematic part of the process is soft inter- 
actions at large distances, where nonperturbative QCD 
comes into the game and manifests itself as the confine- 
ment of quarks and gluons in hadrons and the related 
hadronization process giving the observable hadronic fi- 
nal states in high energy collisions. 

Diffractive processes are sensitive to the details of non- 
perturbative QCD dynamics and provide a way to probe 
the soft and semihard regimes directly. Diffractive events 
are characterized by a leading "target" particle, carrying 
most of the beam momentum, and a well separated pro- 
duced hadronic system. The "gap" in between is con- 
nected to the soft part of the event and therefore to non- 
perturbative effects at a long space-time scale. Diffrac- 
tive deep inelastic scattering (DDIS) offers a particu- 
larly good opportunity to explore the interplay between 
hard and soft physics due to the precise data from the 
electron-proton collider HERA 0, Q • 

DDIS in lepton-proton collisions involves hard scatter- 
ing events where, in spite of the large momentum trans- 
fer Q 2 from the electron, the proton emerges essentially 
unscathed with small transverse momentum, keeping al- 
most all of its original longitudinal beam momentum (for 
reviews on DDIS, see e.g. Refs. [4-6]). The leading pro- 
ton is well separated in momentum space, or rapidity 
y = 1/2 ln(E + p z )/(E — p z ), from the central hadronic 
system produced from the exchanged virtual photon's in- 
teraction with the proton. Thus, this new class of events 



is characterized by a large rapidity gap (LRG) void of 
final state particles. 

Rapidity gaps in deep inelastic scattering (DIS) were 
discovered by the ZEUS and HI experiments at HERA 
3] , but the first discovery of hard diffraction was in pp 
collisions by the UA8 experiment p}. These processes 
had actually been predicted [8[ by combining Regge phe- 
nomenology for soft processes in strong interactions via 
pomeron exchange, with hard processes based on pertur- 
bative QCD. By parametrizing the parton content of an 
exchanged pomeron (or alternatively diffractive parton 
density functions) it is possible to describe the HERA 
data. However, the extracted parton densities are not 
universal, since when used to calculate diffractive hard 
scattering processes in pp collisions at the Tevatron one 
obtains cross sections an order of magnitude larger than 
observed. 

An alternative dynamical interpretation of hard 
diffraction was proposed in Refs. [9[. This Soft Color In- 
teraction (SCI) model is based on the simple assumption 
of soft gluon exchanges leading to color rearrangements 
between the final state partons. Variations in the topol- 
ogy of the confining color fields lead to different hadronic 
final states. 

The SCI model is implemented in Monte Carlo event 
generators, e.g. Lepto for DIS [10]. The hard part of 
the process shown in Fig. Q] is then calculated in the 
framework of perturbative QCD with DGLAP evolution 
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FIG. 1: Schematic illustration of the mechanism for diffractive 
deep inelastic scattering considered in this paper, with soft 
gluon exchanges in the final state. 



of the parton showers in the same way as in inclusive DIS. 
The large momentum transfer means that the hard sub- 
process occurs on a spacetime scale much smaller than 
the bound state proton and is thus "embedded" in the 
proton. The emerging hard-scattered partons propagate 
through the proton's color field and may interact with it. 
Soft exchanges will dominate due to the large coupling 
and the lack of suppression from hard gluon propaga- 
tors. Therefore, the momenta of the hard partons are 
essentially undisturbed — the soft, long distance inter- 
actions do not affect the hard, short distance process, 
and the momentum transfer of the soft exchanges can be 
neglected. However, the exchange of color changes the 
color charges of the emerging partons such that the con- 
fining string-like field between them will have a different 
topology, resulting in a different distribution of final state 
hadrons produced from the string hadronization [11]. In 
particular, a region in rapidity without a string will result 
in an absence of hadrons there, i.e. a rapidity gap. 

The only parameter of this model is the probability 
for a soft exchange, accounting for the unknown nonper- 
turbative dynamics. Remarkably, the SCI model is phe- 
nomenologically very successful in describing many dif- 
ferent processes, both diffractive and nondiffractive [La ]. 
with only a single parameter P ~ 0.5 for this probability. 
Thus, the SCI model captures the essential dynamics of 
diffraction, but lacks a solid theoretical basis. 

To understand better what we can learn from the phe- 
nomenology of the SCI model, we present in this paper 
a detailed QCD-based mechanism for soft gluon rescat- 
tering of final state partons, as illustrated in Fig. [TJ This 
mechanism leads to effective color singlet exchange and 
thereby to diffractive scattering. Inspired by the SCI 
model, the model presented here may be seen as an ex- 
plicit realization of the earlier attempt |13| to understand 
soft gluon exchange in terms of QCD rescattering. Our 
model was initially introduced in a recent letter [14| , and 
is here presented in detail. 

The paper is organized as follows. In Section II we 
briefly discuss the framework of the dipole approach and 
motivate our study. In Section III we consider the kine- 
matics of diffractive DIS. Section IV treats the formalism 
for generalized unintegrated gluon distribution functions 
in the diffractive limit. The explicit calculation of the 
qq dipole contribution to the diffractive cross section and 
analytic approximations used are presented in Section V. 
In Section VI we study the contribution of the qqg final 
state. Numerical results and comparisons with HERA 
data on the diffractive cross section are given in Section 
VII. Finally, in Section VIII we present some concluding 
remarks and an outlook. 



II. DIPOLE APPROACH 

Typical contributions to the diffractive DIS process are 
represented by the diagrams in Fig. [2J In terms of the 
four-momenta q of the photon, and P and P' of the ini- 
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FIG. 2: Typical diagrams contributing to diffractive DIS with 
the leading order qq dipole scattering (left) and the gluonic 
qqg contribution (right) contributing significantly at j3 — > 
(large Mi). 



tial and final proton, the kinematics of the -f*P —} XP' 
process is defined by the variables 



Xb 



Q 2 



Q 2 + W 2 ' 
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Ml' 



x P 
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where Q 2 — —q 2 . The invariant mass of the produced 
system Mx, and the total energy in the -f*P center-of- 
mass system W are given by 



M 



1 ' ] "\ W 2 = (P + q) 2 = ^(l-X B ). (2) 
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The DDIS cross section in general is represented as a 
function of /?, xp, Q 2 and the momentum transfer along 
the proton line t = (P' — P) 2 . Note, that we are working 
in the forward limit of small |i| <C Q 2 , M|. 

Let us consider first the simplest case of the qq con- 
tribution, which is the leading one for small Mx (or, 
equivalently, j3 — >• 1). To compute the diffractive DIS 
amplitude, it is convenient to consider the process in the 
dipole frame [15j, where the deeply virtual photon with 
large virtuality Q 2 and polarization A first splits into a 
quark q and an antiquark q with mass m q , spins a and /3, 
and flavor /, and then the qq dipole with transverse size 
r interacts with the target proton at impact parameter b 
and dissociates into a final state X of invariant mass Mx 
as shown in Fig. [3] The photon splitting into the dipole 
is a QED process and is described by the wavefunctions 
ip^ {z,r; Q 2 ) in the impact parameter space (la. Il6j. 
where z is the fraction of the longitudinal momentum car- 
ried by the quark. The amplitude of the dipole-nucleon 
scattering (denoted as T in Fig. [2]) is the only unknown 
non-perturbative object (for a review, see Ref. |5J). The 
dipole picture naturally incorporates the description of 
both inclusive and diffractive events into a common the- 
oretical framework [lfj, LL7J , as the same dipole scattering 
amplitudes enter into inclusive and diffractive cross sec- 
tions. 

Final states with gluons are suppressed by powers of 
a s . However, if fi becomes small or Q 2 large, the qq 
dipole emits soft or collinear gluons accompanied by 
large logarithms ln(l//3) or InQ 2 which compensate the 



suppression in a s [18J. The qq pair can emit soft glu- 
ons, leading to the dressing up of the quarks, which 
is parametrized by a scale-dependent constituent quark 
mass m|: ff (/i 2 ). In general, the more gluons in the final 
state, the larger the invariant mass produced. The domi- 
nant gluon emission from quarks is described by DGLAP 
evolution [19j and is mostly collinear to the radiating 
quark, so it cannot build up a large Mx- The small 
f3 — >• and large Q 2 — > oo limits can be driven, therefore, 
only by a semihard gluon radiation from the active gluon 
(carrying Xp) giving rise to a gluonic dipole qqg contri- 
bution. These aspects will be discussed in detail below. 

Diffractive DIS at the leading order in a s is described 
by the two gluon exchange diagram shown in Fig. [3l Let 
us first discuss how the longitudinal momentum transfer 
between the qq dipole and the proton can be shared be- 
tween the gluons. The gluon momenta can be Sudakov 
decomposed as 
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-x'P + A', 



x P . 



Cutting the diagram after the first gluon exchange (pick- 
ing up the leading poles only), in the high-energy limit 
xb — > 0, we have 



(q + xP - A ± ) 2 = Afj 



2 
int. 



-2Pq{x B -x)= M; 



2 
int,. 



In the deep inelastic limit Q 2 —> oo, for fixed invariant 
mass of the intermediate system M; n t ~ Mx and Bjorken 
variable xb, we see that M 2 nt <C 2Pq and i~ib, thus 
x' = xp — xp. On the other hand, when Q 2 — > oo and 
M 2 nt fixed, we see that f3 — xp/xp — > 1. So x — > xp 
and x' — > 0, and the first gluon takes all the longitu- 
dinal momentum exchange neutralizing the virtuality of 
the qq system. The latter kincmatical configuration gives 
the leading contribution to the cross section, whereas the 
other configurations with equal momentum sharing be- 
tween the gluons x ~ x' ~ xp/2 are suppressed by extra 
propagators. 

Thus, we consider the asymmetric case with one hard 
(perturbative) gluon carrying most of the longitudinal 
momentum transfer xp, and a number of multiple soft 
screening interactions with total x' <C xp in a color octet 





FIG. 4: Amplitude of the process 7*p 
resummed soft gluon exchange. 



Xp with all-order 



state (in the large N c limit) effectively described by the 
resummed multigluon exchange amplitude as schemati- 
cally sketched in Fig. |4j The first "hard" gluon turns the 
proton into a color octet proton remnant, which then in- 
teracts with a system of soft screening gluons in an octet 
state and finally recombines into a color singlet corre- 
sponding to the leading proton (or system with invariant 
mass close to the proton mass) . These soft gluons cannot 
dynamically resolve a qq dipole of small size r ~ 1//c_l, 
However, we assume that the gluons interact with the 
quark current and not with the dipole as a whole. This 
will be further discussed below in Sec. IV El 



III. KINEMATICS OF DIFFRACTIVE DIS 

Let us first define the kinematics of the process j*p — > 
Xp. Our primary interest is processes with small momen- 
tum transfer t <C Q 2 ■ It is convenient to fix the frame 
of reference in the center-of-mass system of the final 
states, i.e. the outgoing proton with momentum P', and 
the diffractive system X with momentum q' = k\ + k2, 
ki,j_ = -k 2: j_ = kj_. 

The total 7*p c.m.s. energy squared is W 2 = (P' + 
q') 2 = 2P'q' + M x (the proton mass m p is neglected). 
In terms of the cross section variables defined in Eqs. (TTJ) 
and ^ we write 



2P'q' = W 2 



M 2 X 



Q 



,1 — Xp 

xp 



Q 2 

2Pq= — 

xp 



2P'q'. 



Since 2!p < 1, we have W 2 > M x for any j3 and Q 2 , 
so we will first keep M x and then drop it in comparison 
with W 2 whenever appropriate. 

The general Sudakov decompositions of the final 
quark/antiquark momenta fci.2 are 

ki = (1 - z)q + n\P' + fci,_L) k 2 = zq + n 2 P' + k 2: ± ■ 

where n\ = —n 2 . The on-shell conditions for the quark 
and antiquark in the final state 



FIG. 3: Amplitude of the process j*p 
order of perturbation theory. 



Xp at the leading 
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where 



= m 2 + k 2 



is the transverse quark mass 



squared. Applying Eq. @, we get the standard relation 
for the quark transverse momentum 



z{\ 



)M 
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(4) 



and we consider the light quark mass limit m q <C M x . 
The leading contribution to diffractive DIS at HERA 
comes from light quarks, and from now on we do not 
distinguish between their masses and account for them 
by one single mass parameter m e . 

In the diffractive limit Q 2 , M x <§; W 2 , when z and 1 — 
z are not very asymmetric, one has with good accuracy: 
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Let us now define the quark propagators in the pho- 
ton splitting wave function. Due to the condition xp 3> 
^2,3,... the soft screening gluons cannot change the lon- 
gitudinal momenta significantly, but only the transverse 
momenta. Thus, to calculate the hard part of the ampli- 
tude let us first neglect these extra screening gluons. We 
will show below that adding the extra soft gluons leads 
only to phase shifts (and their derivatives) in the trans- 
verse coordinates, which are going be resummed to all 
orders in a s . 

In the chosen frame, the momenta of the exchanged 
hard and the sum of the soft gluons are 



qi ~-x P P' + A ± , q 2 ~A' x , S = \A' ± 



t. 



Let us first attach the hard gluon to the lower quark 
line k 2 ■ Then the denominator of the quark propagator 
between the photon and gluon vertices is 



(fc 2 + <?i) 2 -ml 



-z 2 M 2 x - zQ 2 ^j^ - (k ± + A ± ) 2 - m 
e 2 + (1 - z){k± + A ± ) 2 + zk\ 



(&i 



l-z 

^) 2 
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(5) 



where e 2 = z(\ — z)Q 2 + m 2 The first approximation is 
obtained by substituting M x from Eq. (|3]) in the limit 
xp -c /3 and the second by using the limit A^ -c k±, 
realized when Mx 3> m q ,A±, giving a result valid at 
z< 1. 



When the gluon is attached to the upper gluon line, 
using momentum conservation q = q' + qi + q 2 , we get 
analogously 



(q - k 2 ) - m 



(*U 



A ± ) 2 



(6) 



which is strictly valid at 1 — z <C 1 . It is equal to Eq. ([5]) 
with the exchanges z •<-> (1 — z) and k±_ <-> — fcj_. 

The expressions ([5]) and © will be used for all values 
of z, as is common practice [5]. This is justified in our 
asymmetric case xp ^$> x' because the dominating con- 
tribution to the amplitude comes from the configuration 
that either the quark or the antiquark from the photon 
is essentially on-shell, and the other carries the negative 
virtuality of the photon and then absorbs the hard gluon 
with momentum fraction xp to become essentially on- 
shcll. 

In the limit considered, A^ <C k±, the quark virtuality 
k 2 is conventionally utilized as the factorization scale [ip 
of the process, and is expressed in terms of the energy e 
and the transverse momentum k± as 



H 2 F = e 2 + k 2 L = z(\ 



z)(M x 



+ Q 2 ) 



(7) 



Thus, the hard scale depends on both Q 2 from the space- 
like photon and M x from the time- like final state X. 
Since these can have any values, the QCD factorization 
is complicated and the physics may be different in the 
three cases M x < Q, M x ~ Q and M x > Q- The 
quark propagator (with the hard scale /ip in Eq. (|7|)) 
is antisymmetric, k 2 O —k 2 , with respect to reflection 
between the space-like and time-like regimes, i.e., with 
respect to the interchange M x ■<-> — Q 2 (or e 2 4-> — k 2 ) 
as it should be. 

The next step is to compute the bilinear spinor com- 
binations u{k2,X q )^{X 1 )v{k\,Xq) in the photon splitting 
7* — > qq. The photon polarization vectors in the XP' 
c.m.s. frame have the following general form 

e^(A 7 = ±1) = ^(0,1,±»,0), 
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e^(A 7 =0) 
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(W 2 + Q 2 ,0, 0, W 2 -Q 2 



In particular, by straightforward calculation for the lon- 
gitudinally (L) polarized photon we simply get the fol- 
lowing expression 



u ± (fc 2 )^(A 7 = 0)v^(ki) =* iy/z(l-z)Q . 



(8) 



which is not dependent on the transverse momenta of the 
initial quark and antiquark. 

The transversely polarized case requires a separate dis- 
cussion. Within the dipole picture the diffractive DIS 
process can be basically decomposed into time-ordered 
stages. First, the space-like photon with q 2 = —Q 2 fluc- 
tuates into a qq pair, which is then scattered off the tar- 
get through hard gluon exchange, making the qq-system 
time-like, and finally the on-shell quark and antiquark 



scatter softly off the color background field in the proton 
resulting in a color singlet X-system with invariant mass 
Mx ■ Initially, at the moment of the photon fluctuation, 
the only hard scale is Q 2 , and the transverse momentum 
of a produced quark is expressed through this scale as 
k' ± ~ \/ z (l Z )Q, which is different from the transverse 
momentum k± of a quark in the final state defined in 
Eq. (jj). In particular, the difference between k± and k' ± 
can depend on the actual momentum transfer Aj_ for the 
hard gluon and on the sum of the screening gluons, be- 
cause according to Eq. (£[]) the small variation in k± due 
to the attached Aj_ brings a significant change in Mx 
for small z. The relative coefficient between k± and k'j_ is 
Q/Mx , and has to be taken into account when expressing 
the transversely polarized photon splitting wave function 
through the final state transverse momenta. This physi- 
cal argument agrees well with the kinematics correspond- 
ing to the diagram shown in Fig. 2] for Q ~ Mx ■ In the 
opposite limits Q 3> Mx or Q <C Mx the emission of ex- 
tra gluons significantly complicates the kinematics, and 
this will be considered in detail below. 

As a final result, for the transversely (T) polarized pho- 
ton, we get in the chiral limit 



u±{k 2 )i{\ 1 = ±lK(fci) 



_Q_ 

M x 



2z 
1 -z 



(*?±i*y).(9) 



The spinor signs ± here stand for quark/ antiquark chi- 
ralities, which coincide with the helicity for a quark and 
have the opposite sign of the helicity for an antiquark. 
We do not take into account the "++" and " " com- 
ponents since they are small in the Bjorken limit and for 
relatively light quark and antiquark. 



IV. GENERALIZED UNINTEGRATED GLUON 
DISTRIBUTIONS 



Before calculating the amplitude for the hard and soft 
gluon exchanges, we note that the exchanged gluons all 
originate from the proton color field and should therefore 
be treated through a common description of a general 
gluon density. The first, hard gluon carries the dynamics 
through the longitudinal momentum xpP, whereas the 
soft rescattering gluons carry small momenta x\ -C xp 
but may transfer a color octet charge that screens the 
color of the first gluon resulting in an overall color singlet 
exchange. In this sense, the sum of all soft exchanges acts 
as a single effective gluon exchange between the qq dipolc 
and the proton remnant (see Fig.[5|). 

As an appropriate description of the density of hard 
and soft effective gluons, we use the framework of 
generalized off-diagonal unintegrated gluon distribution 
functions (UGDF), which naturally appear in the k±- 
factorization approach [20 ]. Within this framework 
the coupling to a quark is replaced by an off-diagonal 
UGDF J^ B (x P ,x',A 2 ± ,A' ± 2 ,fi%), absorbing a factor 
Cpa s ([j, 2 olt )/ir and by convention also a gluon propaga- 
tor AT 2 in order to keep the UGDF regular as A 2 ^ — > 0. 
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FIG. 5: The sum of the multiple soft gluon exchanges can 
in an effective description be seen as a single gluon exchange 
between the qq dipole and the color octet proton remnant, 
resulting in the color singlet X-system and a proton in the 
final state. 



The absorbed coupling a s (fi 2 o{t ) corresponds to the cou- 
pling of a screening gluon with virtuality /i 2 oft ~ A^ to 
a quark in the proton, whereas the coupling of the hard 
gluon to the qq dipole is purely perturbative and occurs 
at the hard scale /if- 

Generalized partem distributions (GPD) are not very 
constrained by data. We use a prescription for the gen- 
eralized UGDF [2l| . which works well in the description 
of CDF data on central exclusive charmonium produc- 
tion 22]. This prescription is motivated by positivity 
constraints for the collinear GPDs [23] and can be con- 
sidered as a saturation of the Cauchy-Schwarz inequality 
for the density matrix [241 ]. It incorporates the depen- 
dence on the longitudinal momentum fraction x' and the 
transverse momenta of the soft gluons in an explicitly 
symmetric way, 



Tf ~ y/? g (x P , Ai, &)?„{*, A' ± 2 , Ms 2 oft ) , (10) 

with A± ~ A^_. Here T g is the normal diagonal UGDF, 
which depends on the gluon virtuality A^ and reduces 
to the well-known collinear gluon PDF g(x,/i 2 ) when in- 
tegrated over this virtuality. 

This model (TTUJ) is a factorization of the generalized 
UGDF into a hard part depending on a hard scale fip 
and on xp , thus describing the hard gluon coupling to the 
proton, and a soft part defined at some soft scale /x so f t 
and small x' <C xp corresponding to a number of soft 
gluon couplings. As we will see below, together with the 
factorization in transverse momentum space, the model 
(fT0|) provides a QCD factorization of the diffractive am- 
plitude in momentum space. 

The UGDF F g {x,A 2 Ll fi 2 ) introduced above depends 
on the gluon virtuality, and this dependence is not 
theoretically well-known for small virtualities. The 
UGDF is here modeled using the collinear gluon den- 
sity xpg(xp, fi F ), fixed at the QCD factorization scale 
fip, together with a simple Gaussian Ansatz for the de- 
pendence on the gluon virtuality A]_ as 



x P Tf ~ ^x P g{x P , n%) x'g{x', fi 2 soft ) fa(A 2 ± ), 

f G (A\) = l/(2np 2 ) exp (-A 2 ± /2p 2 ) , (11) 



where the Gaussian width po is the soft hadronic scale. 
As demonstrated below, this scale corresponds to the 
transverse proton size r p ~ 1/po- The Gaussian smearing 
is then interpreted as the result of many soft interactions 
in the bound state proton. The factor y/xp in Eq. (fTTj) 
is absorbed from the hard subprocess part describing the 
coupling of a hard gluon to the qq dipole, and gives us a 
hint th at the off -diagonal UGDF should be proportional 
to ~ y/x P g(xp). 

It is known that at some soft scale /i S oft ~ Aqcd 
collinear PDFs like GRV [25| saturate at small i' < ip, 
so one can introduce a function R g (x', /Jg ft), which is as- 
sumed to be slowly dependent on x' in the case x' <C xp: 



Rgix'^l^xpgixp^DMAl). (12) 

This factor R g (x') contains all the soft physics related 
with soft gluon couplings to the proton. It is interpreted 
as the square root of the gluon density at very small x' <C 
xp and soft scale /ufofr This is a non-perturbative object, 
which contributes to the overall normalization and can 
be determined from data. As will be seen below, the 
prescription (fT2"j) is consistent with the HERA data for 
all available M\ and Q 2 . 

There is a debate about what the power of the gluon 
density in the cross section should be (see [13[ and ref- 
erences therein). When squaring the amplitude contain- 
ing (|12p. this model leads to a linear dependence of the 
diffractive cross section on the gluon PDF. This linear 
dependence is the same as in the SCI model, where it 
describes both diffractive and non-diffractive events, and 
provides a continuous transition between the two types 
of events. 

This is in contrast to the quadratic dependence on the 
gluon density often encountered in two-gluon exchange 
calculations of DDIS @- This arises from another pre- 
scription for the off-diagonal UGDF in the asymmetric 
limit x 1 -C xp and /j. % ^ > Aj_, which in terms of the 
diagonal UGDF reads [26J | 



ff(x P ,x' 



Ai,Ai,/4) 



R g (x')f g (x P ,A 2 



2 \ 

(13) 
where the skewedness parameter 1 R g ~ 1.2 — 1.3 is 
roughly constant at HERA energies, and gives only a 
small contribution to an overall normalization uncer- 
tainty. The factor R g can be approximately taken into 
account in this case by rescaling the xp argument in the 
diagonal UGDF as (23 

ff(xp,x'.Al,Al,p 2 F )^f g (0Alx P ,Al^ 2 F ). (14) 

Using the same Gaussian Ansatz for the intrinsic trans- 
verse momentum dependence as in Eq. (jlip . we see that 
prescription (Ti"3"| leads to a quadratic dependence of the 
diffractive structure function on the gluon PDF. 



The unintegrated gluon density in the form (1101) re- 
duces to the diagonal form (jT5|) in the kinematical do- 
main where x' ~ xp and the hard gluon is soft, fiF ~ 
jUsoft- In this limit there is no QCD factorization, so 
the hard and soft gluons must be taken into account to- 
gether on equal footing. This may be the case at low 
Q 2 and Mjf, when a larger contribution to the cross sec- 
tion comes from relatively soft scales fip < 1 GeV, and 
Eq. (JTT)) reduces to 



■x-p-ff ~ 0.5x P g(0.5x P ,n F )f G (Ai). 



(15) 



where the factor 0.5 appears from the equal momentum 
sharing between active and screening gluons, since xp is 
the sum of all gluon momentum fractions. In this sense, 
the prescription (|10p is more general since it describes 
both x' <C xp and x' ~ xp regimes, and contains pre- 
scription ([TBI as a limiting case. 

Eq. (|15p also leads to a cross section with the gluon 
density in the second power. It is similar to the u R g - 
prescription" (|13p (if the Gaussian smearing like in 
Eq. (TiTl) is adopted), and close to its phenomenological 
form with rescaled argument (I14p . This is more conven- 
tional in the description of the exclusive processes [5| , but 
is valid only for the symmetric case where the two gluon 
exchanges carry longitudinal momentum fractions close 
to each other, x' ~ xp, and are connected to the same 
factorization scale (j,p. This case can also correspond 
to the "no-soft-exchange" approximation, when the soft 
rescattering of the on-shell partons in the final state is 
not taken into account (then R g — 1). 

The above discussion leads to the conclusion that 
whether it is appropriate to use prescription (|12[) or (fT5|) 
depends on what kinematical regime one considers. The 
problem is that x' is not strictly constrained by kinemat- 
ics. Its order of magnitude should be 2 



A 2 



z{\~ z)W 2 



u 2 
Mf 



Xp 



(16) 



Thus the x' <C xp regime is realized in the perturbativc 
limit of large factorization scale /if- In the limit (If — > 
jLtsofti one instead has x' — ¥ Xp, as naturally required by 
matching the prescriptions of Eq. ([12")) and (|15[) . 

To summarize this section, we have formulated a 
framework for the gluon density needed as input to the 
calculation of the diffractive cross section. Since this de- 
scribes soft QCD dynamics in the proton, there are neces- 
sarily some uncertainties. The precise HERA data on the 
diffractive cross section are, however, directly sensitive to 
this gluon density and may, therefore, be used to obtain 
new information about the gluon PDF at extremely small 
x' values and at different scales. 



1 Our function R g (x') is an analog of this skewedness factor. 



2 This is similar to central exclusive production in pp collisions, 
where the screening gluons couple to the triplet/antitriplet 
charges of the proton remnants, which have predominantly equal 
momentum sharing z ~ 1/2, and thus we have x' ~ A3 js. 



V. LEADING-ORDER QUARK DIPOLE 
CONTRIBUTION 

We are now fully equipped to derive the amplitude for 
the dipole-proton interaction. 



A. Hard-soft factorization 

The total amplitude for 7*p ~ > Xp is decomposed into 
longitudinal (L) and transverse (T) parts depending on 
the photon polarization A 7 = 0, ±1, and each part can be 
written as a convolution of the hard and soft subprocess 
amplitudes based on loop integration and cutting rules. 

Starting from the general Sudakov decomposition of 
the total screening gluon momentum q 2 — aq' + bP' + Aj_ 
with a, b <C xp, we can write down the amplitude, for 
example, for transversely polarized photon as 



M- 



A T =± 



d 2 A A 



dadb 



X 5(iki +k 2 + q 2 ) 2 - M 2 x y({P' - q 2 ) 2 ) 
x M hard (A'j_) M soft (AjJ (jfcf ± ik\) . 
The (5-function product can be rewritten as 

which takes care of the integrals over a and b, leading to 

d 2 A 



M 



A 7 =± 



2z 



7i j hard ti ysoft 



±T M x \l 1 - z J (2tt) 2 

x(kf±ikf), 

where, in the frame with q' ± = 0, 

rvi I Zfv-i — — ~~ I run I IKy ) = rv I %t\, 

Analogously, the longitudinal contribution is 

d 2 A 



M. 



A T =0 



i\J z{\ - z)Q 



(2ny 



7i ,rhavd 7i Tsoft 



To calculate the Fourier transform of the total ampli- 
tude we use the convolution theorem 

%) = / f{p)g(q~p)dp -> h(x) = f(x)g(x) 

where / denotes the Fourier transform of /. In our am- 
plitude this convolution is represented by the integral 
over Aj_, while S = y/—i = \A' ± + A±\ plays the role 
of q. Thus, the inverse transformation over the impact 
parameter b is 



M{8) 



d 2 6e- i5b M hard (b) M soft (b) . 



leading to factorization of the amplitude in 6-space as a 
direct consequence of fc^-factorization in impact param- 
eter space. 



B. Hard part 

Consider first the hard gluon coupling to q or q shown 
in Fig. @] in the rotated frame of reference with q± = A± 



and P± 



-Aj_, where kL 



-k' 



k' ± . In what 



follows, we use the relation for quark-gluon vertices in 
the eikonal approximation 

u(k 2+qi )ru(k' 2 )~2k' 2 ^ ?i«*4 

and for the product of the i-matrices in the large N c limit 
we have 

tfjttl — TF$ik5ji • 

The hard part, describing the two possible couplings of 
the hard gluon to the qq pair, can then be written as 



irk^g-ibA^ 



M^ T i (A ± ,k' ± ) = f d 2 rd 2 bM^ d {b,r)e- 

where k' ± is the transverse momentum of a quark in 
the intermediate state (see Fig. [4]), and the Fourier- 
transformed hard amplitudes are given by 

M^ d =iCa 8 ( f x 2 F )^W 3 z 3 / 2 (l-z) 3 / 2 K (sr)V(b,r), 

l = W 2 z 1 ' 2 {l-zf' 2 



M 



hard 



iCa s (p, 




x eK\ (er) - 



f3 yjxp 
V(b,r) ; 



where C = 8ire qy /ira em /N 2 , and _Ko,i are Bessel func- 
tions. The function V(b, r) is the gluon density in impact 
parameter space defined as 



V(b,r) 



1 



u s 0; 
x {e 



soft) 

irAi 



rf 2 A ± 
(2n) 2 

irAj_\ 



^J^ tt (x P ,At 



ibAi 



(18) 



Here, F° n {xp, A z x ) is the generalized UGDF, and fj, 2 oft 
the typical soft scale of the process given by the gluon 
virtuality ~ A\. The factor containing a s (/i 2 oft ) is in- 
troduced in the normalization of the soft part in order to 
compensate its absorption into the UGDF (see above). 
Inserting F° s from Eq. (fT!?)) , we finally get 



V(b,r) 



1 Rg(x') 

a s (p 2 oit ) (2n) 2 



Ib-rl" 



x P g(xp,p 2 F ) 



'-\b+r\- 



(19) 



In the Fourier transformation we assumed slow evolu- 
tion of the QCD coupling ct s (fi 2 o{t ), as is the case in 
the analytic perturbation theory discussed in next sec- 
tion or often assumed by freezing the coupling at very 
small p 2 oit - Thus, in the Gaussian model (fTTj). the un- 
integrated gluon density in the impact parameter space 
V(b,r) is factorized into a collinear gluon density multi- 
plied by an (r, 6)-dependent normal distribution. 
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C. Soft part 

We now turn to the soft subprocess amplitude, which 
can be calculated order-by-order as follows. The soft- 
ness of the color-screening gluons with x\ <C ip implies 
that all intermediate particles are on-shell, and that the 
dipole size r is not changed during the soft interactions. 
Cutting the intermediate propagators we have only the 
phase shifts with the same origin as in Eq. ([T7]). and a 
dependence on the soft momentum exchanges AJ j_ . 

In particular, for one and two soft gluon exchanges (in 
the large iV c limit) we obtain 



-irk' ±M noit =Ae - 



ir\t± 



1 



A'i 



-irA, 



1 



-irk'. 



/-soft 



A 2 



MSOIt 

2! 

-irA', 



^ — irk i 



d 2 A' 2± 



(2n) 2 A'? ± A'L. 



irA- 



-irA', 



1 



where A' 1± = A^ 



A^ and A = 2-ki Cf a s ([i 2 o{t ) with 
Cf — TpN c in the large N c limit. For example, the NLO 
gluonic contribution to the soft part M| oft is represented 
by the four diagrams shown in Fig. [6] 

Fourier transformation with respect to A', leads to 



-irk', jiVsoft 



AI? 



irk I 



i'k' 



M 2 sott = e 



AW(b,r), 
_ irkj ^l 2 W(b,r) 2 



(20) 



2! 



where 



«<>"> = i'»V 



(21) 



Continuing this procedure we see that summing over the 
number of soft gluons in the final state leads to exponen- 
tiation in impact parameter space, so that for the total 
soft subprocess amplitude we finally get 



longitudinal x 1 and transverse momentum transfer A^_ 
at some soft scale ^ so ft- The strong coupling a s (^ 2 oft ) is 
not small in this case. There are several approaches for 
dealing with the Landau pole at low momentum transfer 
(see e.g. Ref. [29j and references therein). We use the 
infrared-finite Analytic Perturbation Theory (APT) [30 ] 
approach to parametrize a s (/i 2 oft ) at /^ 2 oft ~ Aq CD . The 
analytic strong coupling af PT (/z 2 ) is stable with respect 
to the choice of the QCD renormalization scheme, higher- 
order radiative corrections and variations in Aqcd ■ APT 
has also proved to give a quantitative description of light 
quarkonium spectra within the Bethe-Salpeter approach 
31] and DIS spin sum rules at low Q 2 [32]. 

In the one-loop case, the APT Euclidean function Ai, 
i.e. the analyticized first power of the coupling a s in the 
Euclidean domain, is 1301 



,(D 



1 



A 2 



A 2 



Ai 



L = ln 



, A QCD 



AY'(Ai) = - 
Po 

' (23) 

where /3o is the first coefficient of the QCD /3-function. 

Since a significant contribution comes from the phase 
space region with strongly uneven longitudinal momen- 
tum distribution between the quark and the antiquark, 
and where k± is not very large, the diffractive structure 
function becomes sensitive to the model of the strong 
coupling used to calculate a s ° ft , and hence to the typi- 
cal soft scale /x sort of the process. To avoid this problem 
in practice, in the soft regime we are considering, where 
Aj_ ~ Aqcd, we do not extract the value of /x so ft but 
rather fix the coupling at a s ° ft = -Ai(Aqcd) — 0-7. 

Having all the L and T amplitudes, we can now finally 
rotate our frame of reference in the transverse plane to 
one with q' ± — and P' ± = 0, where k± is defined as 
in Eq. (|4]). In this frame, the variables b and r are the 
impact parameter and the qq dipole size in the final state, 
respectively. 



- irk '-LM soft (b, 



-irk I 



(1 



AW(b,T)\ 



(22) 



D. Diffractive structure function 



A similar expression was previously derived in the case 
of scalar Abelian gauge theory in Ref. [28|. Note that 
M soft (b,r) is independent of the photon polarization in 
the soft limit of small A 



i,±- 



As mentioned before, the soft gluon exchanges between 
the final state partons occur at non-perturbatively small 



;*'1,± 



r^ 
















O 


' c 
































o 
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Let us turn to the cross section. We are interested in 
the case when the proton remnant forms a particle in the 
final state with invariant mass close to the proton mass, 
so we have a three particle phase space. We may write 
in general 

, d 3 a*>" p ^ Xp ' 1 



%- 



d 3 q' 2X 1 / 2 (W 2 ,m 2 p ,-Q 

d 3 ki d 3 k 2 d 3 P 



7 ^ (8) (q'-ki 



(2tt) 3 2k° (2tt) 3 2fc° (2tt) 3 2P^ 
x (2tt) V 4 )(g + p-k 1 -k 2 -P<) J2 Wx q x 



FIG. 6: Illustration of resummation of the gluonic contribu- 
tions to the (NLO) soft part M|° ' (two gluon exchanges) in 
the large ^ c limit. 



In the large W limit the flux factor is 
2A 1 / 2 (IF 2 ,m 2 -g 2 ) ~ 2W 2 . The left hand side 



can be transformed to 



E. Physical interpretation and simplification 



, d 3 a^ p ^ Xp ' W 2 d 2 a^ p ^ Xp W 2 f3 2 d 2 cr" f 'P^ Xp ' 



%- 



d 3 q' 



dM%dt 



d/3dt 



The diffractive structure functions have simple relations 
to the corresponding differential cross sections 



?r>W, 



x P F^'((3,x P ,Q 2 ,t) = -^ 
x p F?P(I3,xp,Q 2 ): 



Q*0 d 2 aty Xv ' 



dfidt 

' L,T 



o f 



2 B 1 d 2 af p ^ Xp ' 



47r 2 a em Br 



dfidt 



t=0 



assuming an exponential ^-dependence ~ exp(Bp>t) of 
the cross section on the diffractive slope Bo- The 8- 
functions remove the integrals over P' and one of the 
quark momenta, say, k! , and we get 



d 3 k 2 



z(l-z) 



d 2 al *p^Xp' 


Q 2 n 1 


dfidt 
S(q + P Q -k° 1 - 


4/3 2 (2tt) 5 W 6 



The remaining (^-function removes the integral over k2 Z — 
zW/2. The last phase space integration is rather trivial, 



d 2 ki 



*i 



dkj_ = ttM| I dz(l - 2z) 



Finally, in the full phase space we have to take into ac- 
count an extra factor of two due to the symmetry with 
respect to the interchange z ■<-> 1 — z. 

Straightforward calculation leads to the follow- 
ing expressions for the longitudinal and transverse 
fully-unintegrated diffractive structure functions 
F^\x P ,Q 2 ,P,t): 

1 

x P F° {A) =SQ A M 2 X [' dz(l~2z)z 2 (l-z) 2 \J L \ 2 



x P F° {i) =2SQ i 



(24) 

dz(l-2z){(l-z) 2 +z 2 }\J T \ 2 , 

(25) 



where S = ^„ e 1/(^ 7T ' 2 ^c) sums over light quark charges 
e q , and 



2„ J2 U „-iSb 



J L = ia s (^i 2 F ) I d 2 rd 2 be 



e- trk± K (er) 



V(b,r) 



1-e 



AW 



J T = ia s (n 2 F ) / d 2 rd 2 b 



■2„j2v-i<5b -irk 



eKi(er) 



fx ^t ITy 



V(b,r) 



1 



„AW 



(26) 



(27) 



These are the general expressions of the QCD-based soft 
multiple gluon rescattering model. 



We have now derived Eqs. ([24H27]) . which describe 
the diffractive structure function. These have non- 
perturbative soft gluon exchanges as important ingredi- 
ents, and to calculate these exchanges we have had to 
make some model assumptions. Some of these assump- 
tions have already been discussed above: we treat the cou- 
pling to the quarks using the strong coupling obtained in 
APT and the coupling to the proton remnant using the 
function R g . Moreover, we extrapolate perturbation the- 
ory and assume a perturbative propagator for the gluons. 
The infrared logarithmic divergences in these gluon prop- 
agators, which appear at each order in the resummation 
(see Eqs. (1201211) ). disappear when the gluon exchanges 
are resummed to all orders (Eq. (1221 ). 

There is one additional model assumption, as we will 
explain shortly, but let us first discuss a physical argu- 
ment based on effective field theory principles, or equiv- 
alently, on the uncertainty principle: A gluon with mo- 
mentum k has a "resolution power," or minimal scale of 
an object it can resolve, of order 1/fc. Put another way, 
physics should not depend on scales much smaller than 
the resolution scale. 

The hard gluon in our calculation can resolve the qq 
dipole with transverse size r, allowing us to apply per- 
turbation theory to the hard part. On the other hand, 
the soft gluons do not carry significant longitudinal mo- 
mentum fractions, and only small transverse momenta 
A± ^ 5 = y/—t, so their "maximal resolution" scale is 
b ~ 1/(5. This means that the screening gluons cannot 
dynamically resolve the internal structure of a small qq 
dipole with size r <C b. However, in constructing our 
model, we extrapolate perturbative QCD to the non- 
perturbative regime and assume that the soft gluons cou- 
ple individually to the quark and antiquark, since the es- 
sential point of the dynamics here is the color exchange 
and not the momentum transfer. This is in a similar vein 
to using quark currents in hadronic matrix elements, such 
as form factors. The underlying quark and gluon dy- 
namics is still important even at very low scales (see e.g. 
Ref. [29| for a discussion of this) . In this way there is a 
continuous transition between soft perturbative and soft 
non-perturbative gluons. 

As they stand, the integrals in Eqs. (|26I27|) exhibit un- 
physical singularities in the angular integrations. This, 
however, is because of our model assumption, which so 
far does not fully take into account the resolution power 
argument. Since the gluons are soft, physics should not 
depend on the orientation of the qq dipole with respect to 
the impact parameter. This will regulate the unphysical 
singularities in the angular integration in Eqs. (J26I27I) . 
This will also allow us to evaluate the integrals analyti- 
cally, and we will use this in our calculations below. We 
argue that the resulting expressions, Eqs. (134135)) below, 
are the physically correct expressions for Jl and Jt to 
use in Eqs. (|24I25I) . 

The expression (|19l) can be considered as a model 
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for the unintegrated gluon density in impact parameter 
space. In particular, it defines the probability to probe a 
gluon at impact distance b from the proton center with 
momentum Aj_ ~ 1/6 by a hard qq dipole with small 
size r <C b, where the quarks carry the hard momentum 
k± ~ 1/r. The process is considered at a factorization 
scale equal to the quark virtuality p.'p. The gluons can- 
not resolve scales below the dipole size r. Therefore, the 
gluon density cannot depend on the orientation of the qq 
dipole with respect to b, i.e., on the angle between r and 
b. Also, in this Gaussian model there is no physical rea- 
son for an asymmetry of the UGDF with respect to the 
direction of the vector b. Thus we rewrite our expression 
(fT9|) in the following way 



In the longitudinally polarized case the result of the 
Fourier integration over r is the Hankel transformation 
of Ko(er) r with respect to the momenta k^ and k^ + rj. 
We obtain 



a s {p 2 F ) 



J L ~ Rg -yY ' y/ X pg 



Pi 



,2-u „-i8b -±p 2 b 2 



d"he 



x b 



r 2£ 2 E( K -^r)-{e 2 + kl)Ki K -^) 



2e 2 E - 



e(e 2 + kiy 



(32) 



(kx+rjY 



e (e 2 + (k± + r]) 2 Y 



In the small dipole limit r <C 6 this becomes 

1 



(28) 



V(6,r) 



, ,R g {x')-^^^bre-^ b \ (29) 



a.iuU) " '2n 2 



which will be used below to obtain the formula for the 
DDIS amplitudes. 

We can check our formalism by taking the small cou- 
pling limit a s (p 2 olt ) ~ a s(l^F) "^ 1' wnere we can ap- 
proximate 



1-e 



AW 



-iasi^soh^f 



(30) 



In the longitudinally polarized case, the Fourier integrals 
are then reduced to Hankel transforms, leading to 



in terms of the complete elliptic integrals of the first and 
second kind, K(x) and E(x), respectively. In the forward 
limit of small 5 <C fcj_ we expect fcj_ ^$> rj, rj = \rj\. 

There is a similar simplification in momentum space 
for the hard momentum fcj_ ~ 1/r, i.e., \i±r] ~ k\j\, ne- 
glecting the dependence on the direction of the qq trans- 
verse momentum in the isotropic color field of the proton 
remnant. 

Further, we expand the integrand in Eq. (|32[) in £ = 
k\r\ + r/ 2 <C fcj_, and keep only the leading term in £. 
Taking the last Fourier integral gives 



d 2 be-*e-* 62 



k± + - 



Po 



where 



J L ~ 8R g (x') a s (n 2 F )C F y/x P g t 



e 2 - k 2 



(s 2 



k\) 



Thus, in the limit a s <C 1 our model successfully repro- 
duces the standard leading-order two-gluon amplitude 
(l8f and leads to the correct exponential ^-dependence 
of the cross-section ~ exp(_Bc>£) with diffractive slope 
B D = l/p 2 = 6.9 ± 0.2 GeV- 2 known from HERA data 
[2| • This gives po — 380 MeV, close to the value of Aqqd ■ 
Thus, the Gaussian width po physically corresponds to 
the effective transverse size of the proton. 

However, the strong coupling a s (p 2 Q f t ) is not small in 
the case of small momentum transfers < Aqcd > and we 
cannot calculate the integral in Jj, in general form ana- 
lytically. The soft phase AW is not in general small in 
the Fourier transformation, and in evaluating the Fourier 
integrals in Eqs. (I26I27[) we should not impose the W <C 1 
condition, but rather keep the exponent exp(_4W) with 
imaginary A. This produces an extra phase shift in the 
Fourier transform over r, coming from the soft gluon ex- 
ponentiation in the large N c limit. Employing the "max- 
imum resolution" argument introduced above, we can 
write 



e AW(b,r) _ e - 



1VT] 



V 



<*s(l4ott) C F 



b 



(31) 



WO =/•■■ .<---" ; \l\vp eMlJ-JA, (33) 



I is a modified Bessel function, and v — a s (p 2 o{t )CF- 
The second term is an NLO contribution since it is pro- 
portional to the a s in v, and typically in the forward limit 
t <C 1 GeV and in the hard momentum transfer limit it 
is much smaller than the leading term (we do not consider 
large t, where the whole formalism here does not apply). 
However, this term is the only leading term which sur- 
vives in the limit when both k± — >• and |i| ~ 1/6 2 — » 0, 
so we have to take it into account. 

As regards the ^-dependence, the second term in 
Eq. (|3"3"|) decreases at large t, but not as rapidly as the 
first term. To good approximation, the integral over t in 
the cross section can be written as 

/ dtU{tf ~ pl(k ± + a v) 2 , ct = 0.73 GeV. 
Jo 



The second term must be taken into account when fcj^ is 
< 1 GeV. 

Straightforward calculation leads to the following re- 
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suit for the longitudinal contribution, 

(kl+e 2 ) 



J L ~ Rg^fxpg 



a s {n 2 F )C F U{t) 



2e 3 (kl + e 2 ) 



(34) 



x(8K -^ e z + 7r 2 2 F, -,-;2;-^ e 



1 3 



2'2' 




32e 4 £ 



A- 



and for the transverse contribution, 
a s (n F )C F U(t) 



Jt — Rg^/xpg 



2ekl(kl+e 2 )- 



(35) 



{e 6 + 3k 2 ± s 4 - 2k\e 2 ) E (-^)-(e 2 + k\) 



kl 

r 2 



+ ^ + kl),F 1 (\^2-,- k l)\kl 



(:' + 2Ai.-')A-(-y 



where 2-F1 is the hypergeometric function and -E and if 
are the complete elliptic integrals as above. 



VI. GLUON CONTRIBUTION TO THE 
DIFFRACTIVE STRUCTURE FUNCTION 



In the large- Mx limit, gluon emission may be impor- 
tant. In principle, gluons may be radiated from both the 
qq dipole and the hard gluon. The gluons emitted from 
the quarks are dominantly soft and move collinearly with 
the quarks, and do not significantly change the invariant 
mass of the final system X. Rather, they dress the quarks 
to build up their effective mass m q s , which is, in general, 
a function of the two hard scales Q 2 and M x . This mass 
parameter may be treated as a constituent quark mass. 
In the current work we do not make predictions for m° , 
but instead extract it from data. 

The scale dependence of the effective quark mass in 
processes with two hard scales like the one under con- 
sideration may be complicated. This will be discussed in 
connection with the numerical results in Sec. IVII1 



A. Kinematics 

The small j3 — > limit is, therefore, driven by gluon 
emission from the hard gluon, as illustrated in Fig. [7] 
The kinematics of the process in the XP' c.m.s. frame, 




FIG. 7: Typical diagram for gluon emission in the DDIS final 
state (a), and its asymptotic limit reducing to the "gluonic 
dipole" contribution for fi —t and p± -C k± (b). 



where k 



l,-L 



p±, is given by the Sudakov de- 



compositions 



fci = (1 - z - z')q' + ni P' + k h± , 
k 2 = zq' + n 2 P' + k 2 ,±, 
p=z'q' + n 3 P'+p ± . 



(36) 



with n\ + n 2 + 713 = 0. Analogously to the qq case, 
we obtain the following expression for the invariant mass 
Mx in terms of momentum fractions z, z': 



M 



zk 2 ± + (l-z- z')k 2 ± + (1 - z')m 2 - 2 



x 



z(l- z- z') 



q , P± 

Z' ' 



(37) 



Mx can also be represented in terms of the invariant 
mass of the qq system as 



M 2 - 

M £ = T3i7 + 7(fr7)> (h + k 2 r = M 2 ,. (38) 



p± 



B. Soft gluon resummation and the gluonic dipole 

limit 

The qqg-system scatters off the proton by exchanging 
soft gluons, in the same way as the qq-system above, and 
also here the gluon exchanges can be resummed. We have 
two independent transverse momenta, k[ ± and k' 2 j_, out- 
going from the hard subprocess, corresponding to impact 
parameters t\ and r 2 . Proceeding as for the qq case, we 
obtain for the qqg case 

-irik'i ^ -ir 2 k2 ji j-soft _ „— iriki.j. ^-ir 2 k2.x 
e ' e - lvl l,qqg — e e -T72 



x Ae~ tri ^+Ae 



-ir 2 A.' ± 



■Ac 



where the prefactor A was introduced above for the gluon 
coupling to a qq dipole, and A g corresponds to the case 
with a gluon coupling to a gluon in the qqg-system. By 
explicit calculation of the color factors it can be shown 



that in the large- N c limit A g 



-2 A. This allows us to 



perform the Fourier transformation of the soft part over 
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A'j_ for any number of exchanged gluons and to resum 
them in the same way, as for qq dipole rescattering. For 
one and two gluon exchanges we have 

e -jriki _L e — ir 2 k2 x jy-soft _ e -irik 1) j. e -tr 2 k 3| x 



x A 



W(b,n) + W(b,r 2 ) 



2,999 

A*_ •■ 

2! 



e -iriki ie -ir 2 k 2 



W(b,n) + W(b,r 2 ) 



where W(b, r) is defined above in Eq. (j2"Tj) . Summing 
over the number of soft gluons in the final state leads to 
exponentiation in impact parameter space, i.e., 

e' iri K^ e - ir ^^M aoii {h,T X ,r 2 ) = (39) 

_ e -iriki jie -ir 2 k 2>i ^ — e -A[W(b,n)+W(b,r 2 )] 

Let us now focus on the leading asymptotic behavior of 
the diagram in Fig. EJa) in the limit /3 — ¥ 0. In this limit 
the hard scale of the process fi F ~ Q 2 //3 becomes very 
large. From Eqs. ([37]) or ([35]) we see that the M x -> oo 
limit is realized when z' <C 2 (more precisely z' — >• 0), so 
the invariant mass of the qqg = X system is 



where 



M x 



Mi 



Mi 



A 



» Mi 



z(l-z) 



<M 



x- 



(40) 



(41) 



Consider first the limit where the gluon transverse mo- 
mentum p± is small, such that |ki,_i_| — |1«2._l | 3> |p_l|- 
In impact parameter space this kinematical configura- 
tion corresponds to the diagram shown in Fig. [3(b). 
In this limit the qq pair is very small, i.e., we have 
strong ordering in impact parameter space, which can 
be written as r q <C r g <C b. In color space the qq 
pair can be considered as a single gluon, and we con- 
sider "gluonic dipole" scattering off the target. This is 
consistent with our expression for the corresponding soft 



part (|39p. which in the limit r\ 



i'2 



= r„ reduces to 



- 1 - exp(2^W(b,r g )) = 1 - cxp(-^ g W(b, r 9 )), cor- 
responding to the amplitude for soft gluon-gluon scat- 
tering. This reproduces the conventional gg dipole re- 
sult [18j in the small r g <C b limit, in which the ampli- 
tude of the gluonic dipole scattering differs by a factor of 
N c /Cf — 1/Tp = 2 from the amplitude of the qq scat- 
tering. Indeed, from our model it follows that in this 
limit 



A gg - i 



exp(-.4 g W(b,r s )) 



2ia s ° n C F 



-2 A 



soft. 



(42) 



X 




as compared to Eq. (|3T)|) . 



FIG. 8: Illustration of the qqg contribution to diffractive DIS. 



However, this limiting case cannot give a leading con- 
tribution to the diffractive structure function at large Mx 
because of the smallncss of the transverse momentum of 
the final state gluon p± <C fei,2X- Due to Eq. (|38|) the 
larger gluon p±, the larger invariant mass Mx is pro- 
duced. At the same time, p± can not be significantly 
larger than the quark and antiquark transverse momenta 
&i,2.l- Due to momentum conservation, the maximal Mx 
at fixed z' occurs in the limit p± ~ ki,± 3> k2,±, which 
corresponds to n <C r-i in impact parameter space, lead- 
ing to W(b,ri) < W(b,r 2 ). From Eq. p9]l. this corre- 
sponds to the situation when only the qq component of 
the qqg system scatters off the target with soft part A s °= lt . 
This purely kinematical argument is compatible with an 
observation [33| with respect to models for parton satura- 
tion [34j | , that the qqg and qq dipole contributions should 
saturate to the same value, i.e. A s °^ ~ A s °^ at large 
invariant masses Mx- In particular, this means that the 
scattering of the qqg system off the proton can not be 
reduced to the scattering of the gg dipole. 



Leading qqg contribution to the diffractive 
structure function 



We argued above that the leading qqg contribution to 
the diffractive structure function in the large Mx limit 
comes from on-shell gluon emission from the hard gluon 
as in Fig. [7fa). It is clear from Eq. (|40[) that the rele- 
vant limit Mx — > oo corresponds to essentially on-shell 
gluon emission with z' <C z. The corresponding gluon 
propagator can be only slightly off-shell to give a lead- 
ing contribution to the cross section. In this case the qq 
pair takes most of the longitudinal momentum of the X 
system, and kinematically there is no symmetry with re- 
spect to interchange z' f> 1 - z' in such a qqg system, 
whereas for a qq dipole this symmetry z O 1 — z holds 
explicitly. If one allows the active gluon to couple to the 
qq pair directly, the final state gluon connected to the 
hard quark propagator can not be on-shell, and we get 
an extra suppression of the cross section. Such a "sym- 
metry breaking" in the qqg system does not allow us to 
reduce it to a symmetric gluonic gg dipole and consider 
its soft scattering in the same way as qq scattering. 
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The corresponding physical situation is illustrated in 
Fig. [5] The hard virtual photon first fluctuates into a vir- 
tual qq pair, and the leading configuration is when one 
quark (antiquark) takes most of the photon virtuality 
whereas the other one is almost on-shell. Then the most 
virtual quark (antiquark) emits a (less virtual) gluon, 
which interacts with a slightly virtual sea gluon from the 
proton background field. This last interaction produces 
an essentially on-shell final state gluon, which contributes 
to the final X system. After the first hard gluon exchange 
both quarks have similar and small virtualities and scat- 
ter off the proton background field. 

In order to calculate the qqg contribution to the diffrac- 
tive structure function we include a DGLAP splitting of 
the hard gluon (with longitudinal momentum fraction 
Xp) into two gluons — one carries momentum fraction 
ZgXp and couples to the hard part, and one is on-shell 
and contributes to the final state in 7*p c.m.s. frame as 
shown in Fig. [8j The diffractive structure function cor- 
responding to the qqg contribution can be then written 
as (see e.g. |35| ) 



x P F D{4) 
• hpr qqg 



1 



QXgClZg ~ \^s\^g/ 

*gg\ z g) 



t g + m 2 g 



2- 



xpF q D g i4) (43) 



where the integral is regulated in the infrared by the ef- 
fective gluon mass m g ~ Aqcd in the gluon propagator. 
The factor 1/N 2 is due to averaging over the color indices 
(in the large N c limit) of the extra gluon contributing to 
the color singlet X, and P gg (z g ) is the gluon-gluon split- 
ting function 



P gg {z) = C A 



1 -z 



1 



+ z(l-z) 



(44) 



Since the qq contribution is dominated by transverse pho- 
ton polarization, in our formulation the same is true of 
the qqg contribution. 

One could also include more gluons in the final state 
by applying DGLAP evolution of the gluon density, and 
partially populate the rapidity gap by extra hadronic ac- 
tivity from the hadronization of gluons emitted from the 
hard gluon in the same way as in Monte Carlo simula- 
tions. This would lead to a model describing a smooth 
transition between diffractive and non-diffractive final 
states. 

Let us finally comment on another approach to re- 
summing multi-gluon exchange, which results in similar 
eikonal factors [1 — exp(. . . )] in the amplitudes. This 
approach was developed by Hautmann, Kunszt, and 
Soper [36| (HKS) and by Hautmann and Soper [37] (HS), 
and is applicable in both inclusive and diffractive DIS. 
This approach is similar to ours, employing factorization 
and resummation of soft ^-channel gluons. In diffractive 
DIS, the incoming partonic dipole is assumed to move 
closely together in the transverse plane before interacting 
with the color field of the proton. In the HKS/HS ap- 
proach all exchanged gluons are treated on the same foot- 
ing. These gluons collectively carry color singlet charge 



and are resummed using a Wilson line. In our approach, 
we use conventional /^-factorization in terms of the unin- 
tegrated gluon distribution, and we additionally factorize 
the "hard" gluon, which carries most of the momentum, 
from the rest of the exchanged gluons, which are much 
softer (x' <C Xp) and are resummed to all orders. The 
resummed gluons collectively carry color octet charge, 
which combined with the first gluon is required to form 
an overall color singlet exchange. It would be interesting 
to examine the connections between the two approaches 
further. 



VII. NUMERICAL RESULTS 

The HERA data dH on DDIS are given in the form 
of the reduced cross section 



xpa?W 



x P F D(3) 

Xp qq,T 



2-2y + y 



2 y „D(3) . r D(3) 

2 X P F qq:L + XpFqqg 



expressed in terms of the diffractive structure functions 
F L ^ (xp,Q 2 , j3). The momentum transfer t is inte- 
grated over since in most of the data the leading proton 
is not observed, and diffraction is cquivalently defined 
through a large rapidity gap. The kinematical variable 
y = Q 2 /(sxp) < 1, where y/s = 318 GeV is the center- 
of-mass energy of ep-collisions in HERA. In Fig. |H] we 
compare the latest ZEUS data [2] with the numerical 
evaluation of our model. A generally very good agree- 
ment is found, but this needs to be discussed in detail in 
order to gain understanding of the dynamics involved. 

As discussed in Sec. IIV1 we need the generalized gluon 
distribution function in the proton, and use the prescrip- 
tion in Eq. ([T2"T) for the UGDF. This reduces the problem 
to an input of a standard parametrization of the gluon 
density in the proton, i.e. xg(x, up). Here we mainly use 
the recent CTEQ6L1 parametrization (38[, which is in 
leading order and thereby consistent with our treatment. 
Below we also consider other parametrizations to illus- 
trate the uncertainty at very small x and factorization 
scales \xf. The minimum factorization scale \xf is fixed 
to be Hp mm — 0.2 GeV giving rise to a minimum possi- 
ble fraction of the quark longitudinal momentum z m { n in 
the phase space integral. 

The physical parameters that are fixed, are the "soft" 
coupling a s ° lt = A\ (Aqcd) — 0.7, obtained from 
infrared-finite analytic perturbation theory (see Sec. lVCp 
and used for the coupling of the soft screening gluons, 
and the gluon mass m g ~ /i so ft — Aqcd adopted as the 
infrared regulator in the gluon propagator in Eq. (|43|) . 

Fixing a s s ° lt and m g , the only free parameters in our 
model are to° s and R g , representing different soft effects 
that cannot be calculated or safely estimated. The con- 
stituent quark mass mt , which enters the kinematics in 
Sec. Ill, accounts for the soft gluon radiation from the 
qq dipole and corresponds to forming dressed quarks be- 
fore hadronization. The soft part R g of the off-diagonal 
UGDF in Eq. (fi"2"]) can be identified with the square root 
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FIG. 9: The reduced cross section xpa r ( ' (xp, /?, Q 2 ) as a function of xp for different values of Mx and Q 2 . The latest ZEUS 
data 01, from diffractive deep inelastic scattering events with a large rapidity gap, compared with our model using for the 
gluon density in the proton the CTEQ6L1 (full line) parametrization [3g| and at low x and Q 2 also the GRV94 (dotted line) 
parametrization. 
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of the "soft" collinear PDF defined at some x' <C xp and 
/i so ft which represents the soft scale of the color screening 
gluons. The sensitivity to these parameters is discussed 
in the following. 

The shapes of the model curves are in quite good agree- 
ment with the data, except for a few points at extremely 
small xp < 5 x 10~ 4 , and small scales Q 2 < 5 GeV 2 
and Mx (see the upper right corner of Fig. [SJ. Here, 
we are in the kinematical domain where the uncertain- 
ties in the parametrizations of the gluon density of the 
proton become extremely large, as illustrated in Fig. [TO] 
As can be seen, for x < 10 -3 there are substantial differ- 
ences between the different gluon parametrizations and 
the differences become huge for x ~ 10~ 4 and Q 2 < 1 
GeV 2 . The reason is that there is no data from inclusive 
DIS or other processes that can measure the gluon den- 
sity directly in this domain. This gluon PDF uncertainty 
strongly affects the calculated diffractive structure func- 
tion at small quark fractions z and/or small Q 2 and Mx , 
where [ip may drop below 1 GeV. In this case x' ~ xp, 
due to Eq. (ITC|) . so our basic assumptions and QCD fac- 
torization itself become less reliable. In principle, the 
diffractive DIS data can be utilized for selecting the best 
gluon parametrization among those available in the liter- 
ature or, even better, for making new gluon parametriza- 
tions including the data on DDIS, which depend directly 
on the gluon density. 

Another signature of such uncertainties is the Q 2 be- 
havior of the soft part of the UGDF, i.e. R g (Q 2 ), which 
is shown in Fig. [TT] for different PDFs. As can be 
seen, using either the typical leading order CTEQ6L1 
PDF [38], which decreases at small xp and fip, or the 
more regular but older GRV94HO PDF [25[ (see com- 
parison in Fig. I10[) to perform the fit of our model to 
ZEUS data results in quite different fitted R g (Q 2 ). At 
higher scales, Q 2 > 16 GeV 2 , the soft factor is quite sta- 
ble at R g ~ 1. However, at Q 2 < 5 GeV 2 , the diffrac- 
tive cross section calculated with CTEQ6L1 is underesti- 
mated by almost an order of magnitude, and in order to 
get the correct normalization the fitted R g value grows 
significantly. This is mostly because of the strong sup- 
pression in CTEQ6L1 at small xp < 5 x 10~ 4 at scales 
/if, ~ 1 GeV 2 . In contrast, the fit with GRV94HO, which 
does not decrease at small xp, leads to a more stable 
behavior at low Q 2 , such that R g essentially becomes an 
overall normalization constant close to unity. 

In order to illustrate how uncertainties in the PDFs 
and in the UGDF prescriptions affect the xp-dependence 
in comparison with data, we compare our model to the 
data using both CTEQ6L1 and GRV94HO in the UGDF 
prescription of Eq. (|12[) . This is our "normal" prescrip- 
tion that gives a linear dependence of the cross section on 
the gluon density. For comparison we also use CTEQ6L1 
with the "old R g " -prescription defined in Eq. (IT5|) . which 
makes the cross section depend on the square of the gluon 
density. Fig. [T3] shows the results in the bins of interest 
with small scales Mx and Q 2 . The "old R g " -prescription 
leads to an order of magnitude too small diffractive struc- 



ture function at all Q 2 , which cannot be explained by 
the expected normalization factor R g of order unity in 
Eq. (fl"3"|). The corresponding curves in Fig.[T2"lIrave there- 
fore been normalized in order to compare the xp-slopes. 
These slopes are in reasonable agreement with the data, 
but there is a tendency for a too large curvature gener- 
ated by the squared gluon density, in particular at large 
Q 2 . 

For the curves with linear gluon density, the curves 
with GRV94HO lead to better slopes at the smallest Q 2 
and Mx than those with CTEQ6L1, but at higher scales 
they become too steep and cannot describe data. This is 
not surprising since the old GRV parametrization from 
1994 docs not take into account later data from HERA 
and elsewhere, but it provides an interesting alterna- 
tive due to its more regular behavior at very small x at 
low scales. The curves fitted with the recent CTEQ6L1 
parametrization have better xp-slopes at higher scales 
and this is therefore the main alternative in Fig. [SJ in 
spite of its shortcoming at the very lowest xp points. 

The remaining free parameter to discuss is m° , the ef- 
fective mass of the quark and antiquark in the A-system 
which is used in kinematical relations. In Fig.[T2]we show 
fitted values of m* at different scales Mx and Q 2 . The 
diffractive cross section itself is not very sensitive to m^ ff , 
which therefore only varies within the physically reason- 
able interval Aqcd < ™? g ^5 1-3 GeV. Thus, it is mostly 
of nonperturbative nature and can be interpreted as a 
constituent quark mass. Nevertheless, m° depends on 
both Q 2 and Mx, indicating that both scales contribute 
to generating softer gluon radiation. This dependence is, 
however, non-trivial. 

Indeed, a larger invariant mass Mx provides a larger 
phase space, which may accumulate more soft collinear 
gluons, leading to a larger effective quark mass. On 
the other hand, a harder photon (large Q 2 ) can probe a 
quark at smaller distances, so the Q 2 -dependence of the 
quark mass obeys renormalization group evolution, i.e. 
TO° should decrease at larger Q 2 . These two effects are 
indeed observed in the description of data (see Fig. (fTB")) 
for Mx < 10 GeV. At larger Mx the situation changes 
somewhat due to more hard gluon radiation contribut- 
ing to Mx (in particular, the qqg contribution becomes 
important). 

We now investigate the role of the different contribu- 
tions to xp(j r in Eq. (|45p. i.e. qq from longitudinally 
and transversely polarized photons and the qqg contribu- 
tion. In our results shown in Figs. 191 and IT1Z1 above, they 
are all included. We find, however, that the leading order 
gg-dipole contribution dominates in all bins of Mx and 
Q and is enough to describe all data for j3 > 0.2, be- 
low which the qqg contribution becomes significant and 
can be approximated with its leading part calculated via 
DGLAP splitting of the first, hard gluon. 

Fig. [T4l shows the /3-spectra for the different contribu- 
tions. They all vanish in the limits of j3 — > 0, correspond- 
ing to large Mx, and (3 — > 1, corresponding to small or 
vanishing Mx with production of resonances, which is 
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FIG. 10: The dependence on momentum fraction x for the gluon density in the proton at the scales Q 2 — 1, 2 and 5 GeV 2 
given by the standard parametrizations CTEQ6L1 [H, GRV94HO [Hj], MSTW2008LO, and MSTW2008NLO [3J1. 



not taken into account here, or no available phase space. 
The transverse qq contribution dominates over the other 
contributions. The longitudinal qq contribution is always 
small, although it becomes slightly larger at smaller Q 2 
scales. The gluonic qqg contribution becomes relatively 
larger both at high Q 2 and small (3, where it gives an 
important contribution that must be taken into account. 



VIII. CONCLUSIONS AND OUTLOOK 

We have in this paper developed a proper QCD frame- 
work for diffractive hard scattering, which contains both 
hard and soft dynamics. The hard part produces a well- 
defined state of emerging partons, and the soft part is the 
rescattering of these partons with the color field of the 
proton remnant. We have demonstrated that, by taking 
the Fourier transform from momentum space to impact 
parameter space, the overall amplitude can be factorized 
into separate amplitudes for these hard and soft parts. 
This provides a substantial simplification for the calcula- 
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tion and is consistent with the physical insight that soft, 
long distance processes cannot affect the hard process 
occurring on a short distance scale. 

The hard part is calculated using perturbative QCD, 
in the same way as for inclusive DIS. A perturbative hard 
scale is provided by the photon virtuality Q 2 and invari- 
ant mass Mx of the diffractive system, and the process 
thus occurs at a space-time scale much smaller than the 
proton size. For small x, the hard subprocess "f*g — > qq 
dominates. This process is mediated by a single gluon ex- 
change taking most of the longitudinal momentum trans- 
fer, and leaves a proton remnant consisting of the three 
valence quarks in a color octet state. The proton remnant 
carries most of the beam momentum, and is therefore well 
separated in rapidity from the qq system. 

The soft part of the amplitude accounts for the rescat- 
tering of the qq pair (in a color octet state) with this 
remnant. This rescattering is dominated by multiple ex- 
changes of soft gluons, which have larger couplings and 
less propagator suppression. The result is a negligible 
change of the momenta of the emerging partons, but an 
important change of phase is picked up — this is the 
essence of the eikonal approximation. We find that sum- 
ming over an arbitrary number of exchanged gluons leads 
to exponentiation of the soft amplitude, which can be 
written in a closed analytic form free of infrared diver- 
gences. The color exchange, treated in the large- A^ c ap- 
proximation, leads to an overall color singlet exchange 
between the qq dipole and the proton. These two color 
singlet systems then hadronize independently separated 
by a gap in rapidity as characteristic signature of diffrac- 
tive scattering. 

By invoking physical considerations based on the un- 
certainty principle, which limits the possible resolution 
of small momentum transfers, we obtain simplifications 
of the otherwise complicated angular relations in the im- 
pact parameter space. In essence, the orientation of the 
gg-dipole relative to the proton color field is physically 
irrelevant and can be averaged out. 

In addition to the leading order contribution from the 
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FIG. 12: The reduced cross section xpa r (xp, /3,Q 2 ) as a function of xp in the region of low scales Mx < 3, 5, 8 GeV 
and Q 2 < 12 GeV 2 , where the lowest a^p-values are reached. ZEUS data Q| are compared with the results of the model using 
CTEQ6L1 (solid line) and GRV94HO (dotted line) in the "square root" UGDF prescription JT2), and with CTEQ6L1 in the 
u R g " prescription (|15l) (dashed line, normalization adjusted for easy comparison with the other curves). 



qq-dipole, we have also included the next-to-leading con- 
tribution qqg with an extra gluon in the final state. Here, 
we find that the most important contribution is emis- 
sion of this gluon from the exchanged hard gluon (in 
7*g — > qq), which can be well approximated by leading 
logarithmic DGLAP emission. 

Numerical evaluation of the analytical results gives 
good agreement with the precise HERA data on the 
diffractive deep inelastic cross section. The qq contri- 
bution is indeed dominant, but at /3 < 0.2, the qqg con- 
tribution is important. At very small xp < 5 x 10~ 4 
and scales fip ~ 1 GeV 2 the gluon density in the proton, 
which is used as input in our calculation, is very poorly 
known and gives a complication in the comparison with 
the few HERA data points in this extreme region. Stan- 



dard up-to-date parametrizations have a too low gluon 
density in this x, fip region, whereas, e.g., the old GRV94 
gluon density does better. Since the diffractive cross sec- 
tion depends directly on the gluon density, and not only 
indirectly via DGLAP evolution as for inclusive DIS, one 
here obtains an interesting possibility to constrain the 
gluon density at very small x. 

Having demonstrated that our theoretical formalism 
for DDIS does describe HERA data, one may then ex- 
tract the part describing the multigluon exchange process 
and apply it to other hard scattering processes. This soft 
rescattering description ought to be universal, due to the 
factorization of the hard and soft amplitudes. Thus, one 
may apply it together with hard processes in pp collisions 
at the Tevatron to describe the different hard diffractive 
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FIG. 13: The effective quark mass parameter m e q extracted from the fits of the HERA data shown in Fig. [51 using the 
CTEQ6L1 gluon parametrization in the "square root" prescription (|12[) . In the left panel, as function of Q 2 with curves for 
Mx — 20 GeV (dashed) and Mx = 5 GeV (solid). In the right panel, as function of Mx with curves for Q 2 = 5.5 GeV 2 
(dotted), Q 2 = 16 GeV 2 (dashed) and Q 2 = 110 GeV 2 (solid). 
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FIG. 14: /3-dependence of the different contributions to the diffractive structure function, Eq. (|45[) . at Q 



and Q 2 = 40 GeV 2 (right). Solid lines represent the transverse part xpF c 



D(3) 

qq,T ' 



dashed lines the longitudinal part XpF '' L 



12 GeV 2 (left) 
D(3) and 



dotted lines the gluonic qqg contribution XpF ^ . As usual, the prescription of Eq. (|12|) for the unintegrated gluon density is 



used with the CTEQ6L1 parametrization. (The quark effective mass is here fixed at mj 
these scales.) 



0.9 GeV, being the mean value at 



processed observed there, and then go to the higher ener- 
gies at the LHC. One may also apply it to more detailed 
observables in diffractive DIS, such as diffractive dijets 
or diffractive vector meson production. 

However, not only diffractive processes are of inter- 
est. The Soft Color Interaction (SCI) model discussed 
above has previously been successfully applied to both 
charmonium production and B-meson decays [4y|, and 
one may expect the model presented in this paper to 
have interesting applications in such processes too. More- 
over, the multi-gluon exchange mechanism will also affect 
the underlying event, since it effectuates color exchanges 
that modify the color-string topology and thereby the 
hadronic final state after hadronization. The underlying 
event is important in its own right to understand non- 
perturbative QCD dynamics, and also for understanding 
of inclusive events when subtracting the Standard Model 
background in searches for new phenomena at LHC. 



Finally we note that in deriving the theoretical formal- 
ism presented here, we have not used any assumptions or 
results from the previous Soft Color Interaction (SCI) 
model. Our new formalism stands on its own, based on 
QCD theory and basic physical arguments. The formal- 
ism can, however, explain why the simple SCI model has 
been so successful in describing data on diffractive hard 
scattering and other phenomena. The assumptions of 
the SCI model as well as its major features are essen- 
tially what comes out as results of the present paper. Of 
course, our new formalism has a richer dynamical struc- 
ture and we will therefore attempt to improve the Monte 
Carlo implementation of the SCI model by replacing its 
fixed probability for soft gluon exchanges with a mecha- 
nism based on the above amplitude for the multiple soft 
gluon exchanges. This will introduce a non-trivial depen- 
dence on the kinematical variables, giving a new level of 
event-to-event variations. As usual with full event sim- 
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ulation using Monte Carlo, this will give access to more 
detailed studies of both the employed theoretical model 
and its comparison to data in terms of the indicated more 
elaborate observables. 
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